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Abstract 

This paper gives a PDE for multi-time joint probability of the Airy process, which generalizes Adler 
and van Moerbeke's result on the 2-time case. As an intermediate step, the PDE for the multi-time joint 
probability of the Dyson Brownian motion is also given. 

1 Introduction 

The Airy process can be defined purely stochastically as the limit of the Dyson Brownian motion, as we are 
going to do later. However, it also appears in various statistical physical models, such as the polynuclear 
growth process [H] , [7] and the Domino tiling model [H] . Since the Airy process is stationary with continuous 
sample path [12j . we can pick any time t and consider the probability of all particles being in (— oo, u), denoted 
by P(w), and find that the probability is given by the GUE Tracy- Widom distribution [13] 

P(m) (1) 

with q{s) the solution of the Painleve 11 equation 

for s — > CXD, 



q"is)^sqi.s) + 2q\s), q{s) \ ^ (2) 

I y/—s/2 for s — oo. 

In their study of the joint probability for several times of the Airy process, Prahofer and Spohn [12j posed 
the problem to find a PDE for the joint probability. Adler and van Moerbeke [2] solved the problem for the 
2-time case, and assuming a plausible conjecture of the boundary condition, got the asymptotic expansion 
of the probability function ¥(t,u,v), which is the probability that all particles are in (— oo,m) initially and 
in (—CXD, v) after a large time t. Their solution was obtained by a previous result of them on the spectrum of 
coupled random matrices [!lj. They regarded the joint distribution for the Dyson Brownian motion of 2-time 
as a r function of the two-Toda lattice, and construct a PDE with variables in times and boundary points 
of the Dyson Brownian motion as a consequence of identities for r functions and Virasoro identities specific 
to the situation. Then they got the PDE for the Airy process by taking the limit. 
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This paper generalizes their result to the multi-time case, and the technical heart is the same identity 
for T functions, although in the generalized case we need more elaborate work to fit differential operators in 
times and boundary points of the Dyson Brownian motion into the structure of two-Toda r functions. 

After the description of the problem, we state the PDEs for both the Dyson Brownian motion with finite 
particles and its limit, the Airy process with infinitely many particles, and an example for the 3-time (m = 2) 
case for the Airy process. Section [2] derives the result for the Dyson process and section [3] derives the result 
for the Airy process by taking limit. 

1.1 Description of the model 

The free Brownian motion process is determined by the transition probability distribution 

P{t, X, X) = e W]3- , (3) 

v/(2^ 

where X and X are initial and terminal coordinates of the particle, and P is the diffusion constant. The 
probability distribution P(t^ X , X), as a fimction of t and X, satisfied the diffusion equation 

^.^^P. (4) 

dt 2/3 dX^ ^ ' 

If we add a harmonic potential pX^/2 to the process, then the probability distribution P{t,X,X) satisfies 
(see e.g. [5]) 

and the process is determined by (c — e"''') 

1 (X-cX)^ 

P{t,X,X) = ^==e (l^™. (6) 
V pP 

While the free Brownian motion process is dispersive, the Brownian motion process in the harmonic potential 
well has a stationary distribution 

nx) = ^==. (7) 

Now we can define the Ornstein-Uhlenbeck process |llj of an n x n Hermitian matrix B, in which all 
the n? real variables — n for real diagonal entries, n{n — l)/2 for the real parts of off diagonal entries, and 
the other n{n — l)/2 for the imaginary parts of them — are in independent Brownian motion in harmonic 
potential wells. The p for them are uniformly 1, and (3 are I for the n diagonal variables and 2 for the 
n(n — 1) off diagonal variables. Therefore for i, j in {I, . . . , n}, (c — e^*) 

P,,(t,Bu,Bu) = , e''"'V^"' , 

^7r(l-c2) 

P,si (i, , ) - ^7^=^e , (8) 
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and we can write the joint transition probability distribution as 



2—1 l<i<j<n 



We consider the multi-time transition function with the initial state at to = 0, the terminal state 
B„i and a series of intermediate states Bi, . . . , B^-i, and the time between state Bq and Bi being ti, if we 
denote 

h ^ = (10) 

ti-ti^i i = 2,...,m, 
and 

Q = e-«% (11) 

then 

P{ti,...,tm,BQ,...,B„,) = G-^\{e ' . (12) 

1=1 

The Ornstein-Uhlenbeck process has a stationary distribution 

P(B) = C7-ie-T^^'. (13) 

Since the Ornstein-Uhlenbeck process is invariant under the unitary transiformation, we define the pro- 
cess of the eigenvalues as the Dyson Brownian motion process |4j , whose multi-time transition probability 
distribution is (0 = < < • • ■ j < im) 

The transition probability of the n x n Hermitian ma- 

/ s trix with eigenvalues initially A^"^ — (X\^\ . . . , 
P(ti,...,t„„A(°),...,A(™)) = ^ VI ' ' " (^4) 

and A^ after time ti, > after time t2 ■ ■ .and X^"^' 

after the total time t^- 

2 

If we change the coordinates of the M" space of n x n Hermitian matrices in to the eigenvalue-angle 
coordinates Ai, . . . , A„, 6*1, ... , 0„(„_i), with the Jacobian identity (see e.g. [TD]) 

n n n{n — l) 

lldxu n {d'fi{x^J)d^{x^/j)^ViXfY[dX^ Y[ de„ (15) 

2—1 l<z<j<n 2—1 2—1 

where V{X) — Y[i<i<j<ni'^i~-^j) Vandermonde, we find the explicit formula for P(ti, . . . , t^, A'-'^-' , . . . , A'-™-'): 

p(ti,...,t„,A(°),...,A(™))=i /... /ne"^"' n^(^''^)'n^^^^^ (i^) 

i = l 2=1 2=1 

where S*-*^^ does not appear in the integral since the transition probability is independent of 0'"' for the 
unitary invariant property. 

By the Harish-Chandra-Itzykson-Zuber (HCIZ) formula [3] 

(7(22) V{x)V{y) 
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where X = diag(a:i, . . . , Xn) and Y — diag(?/i, . . . , ?/„) are diagonal matrices, we can evaluate the multi-time 
transition probability density as 

P(ti,...,t™,A(°),...,A(™)) = lnA(°))-V(A(™))ndet ' ' ) 

1=1 ^ ^ 



If we take the initial state with eigenvalue A^°-' from the stationary distribution ([13]), which is 

P(A(")) = i^(A(0))2e-2:r^i^''"', (19) 

We get the multi-time correlation function in the stationary Dyson process 
P(ti, . . . , i™, A("\ . . . , A(™)) =P(A("))F(^i, . . . , t„, aW, . . . , Af")) 



1 ™ / 2o, 



(20) 



! = 1 \ / 

If we want to find the probability of all Af 's being in ~ (aJ'^Oj^) U ••• U (a2n_i, a2r,), with 
—oo < a^'^ < 02 < flg ^ < • • • < 02^!^ < oo, for / = 0, 1, . . . , m and « = 1, . . . , 7i, which is 

■■■ / F(ii,...,t„.,A(°),...,AW)nn'^4^ (21) 

we can simplify it by the symmetry of A^'^* , ■ • ■ , aI'-* , for alH = 0, . . . , m and get 
pDy^°"(t,;af)) = 1 / . . . / y(A(0))F(A(™)) 

i=l i=l i=0 fc=l 

We are going to give a PDE satisfied by log pDyson ^\\]^ variables ti and af . 

The Airy process can be defined as the limit of the Dyson process on the edge [5]. As n ^ oo, we can 
prove that the right-most particle in the Dyson process is almost surely around \/2n with the fluctuation 
scale [H], [2]. If we take the rescaling 

ti ^n^'hu (23) 

A«=V2nV6(A«_V2;i) (24) 

a^P =^T?''^(ai ^ ~ V2^), (25) 

then for fixed U and af\ pDyson 

converges to a function defined by the Fredholm determinant of a matrix 

integral operator [5], 0, [H], [7] 

JimpDy-"! ,^^,^/„v3 =pA-''(ii,...,t„,al°\...,4:2) = det(/-(x-if.^X,^)i<^,<™), (26) 
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where Xi is an indicator function defined as 



f otherwise. 



and (Ai is the Airy function) 



KfAx,y)=i 



Ai(x + z)Ai(?; + z)dz if i = j, 

e-^(*--*^)Ai(a; + z)Ai(?; + z)dz if i > j, (28) 
- e^(*J'^*')Ai(a; + z)Ai{y + z)dz if i < j. 



Then we can define the Airy process, which contains infinitely many particles by the probability function 
(PS)) . Furthermore, we are going to give a PDE satisfied by r with variables U and 

Remark 1. To make the definition meaningful, we need a^''' to be —oo for all I. Otherwise the left hand 
side of ([^5]) is and the right hand side is not well defined. 

1.2 Statement of main results 

With notations defined in subsection (|l.ip . we define differential operators {I = 0,1, . . . ,m) 

2ri „ 2ri „ 

^''^ = E4"-%' (29) 

i=i 9a\ fe^i da] 

if all af'' are finite; otherwise we drop the a^''' (resp. 02^^) part if af^ = —00 (resp. Oj^^ — 00). And then 
denote 

m 

yli =^e~*'i?'^i, (30) 
/=o 

m 

Si =^e*'-*"Z?''i, (31) 
y^2=^e-'*'i?''' + ^(l-e-2*^)^-e-^*™, (32) 

mm o 

82= J2 e2(*'-*")i?''2 _^ ^(g2(t,-t,„) _ e-^*-)^ - 6-2*-. (33) 

/=0 1=1 ' 

We now state 

Theorem 1 (Dyson Brownian motion). Given ti, . . . ,tm, the logarithm of the joint distribution for the 
stationary Dyson Brownian motion pDyson ^gyjj^g^ j,^ (abbreviated as logP„J satisfies a third order 

non-linear PDE in times and boundary points of t/^'^ 

B2A1 logP„ ^ -^agilogP^ 
^6iy^ilogP„ + 2ne-*- ^y^i6ilogP„ + 2ne-*'"" ^ ' 
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Similarly with the notations 



in 

1=0 
m 

1=0 

m 

Vi = = ViL - = - 2ti)D 

1=0 

m 

rn 
m 
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1=0 



1=1 



1=1 



(35) 
(36) 
(37) 
(38) 
(39) 
(40) 
(41) 
(42) 
(43) 
(44) 



we state the result for the Airy process 



Theorem 2 (Airy process). Given ti, . . . ,tm, the logarithm of the joint distribution for the Airy process 
pAiry fi^jijiQfi ( |g61) (abbreviated as logP^ satisfies a third order non- linear PDE^ in times and boundary 
points of t/^'-* 

V^[Ei + I?3 + logP - VVi[E + V2+ Ti] logP - 2VilT>irVi logP = {V^ logP, VVi logPjp. (45) 

In the case of to = 1, our result agrees with that in 2^. Especially, if = {~oo, u), J/'^^-' — (oo, v) and 
denote ti = t, then the result for logP^"'^(i, u, v) is 

Corollary 1 (|2j). The logarithm of the 2-time joint probability for the Airy process P^"^^ (t, u, v) ( abbreviated 
as logP^ satisfies a third order non-linear PDE in variables u, v and t 



{v — u) 



d_ d_ 

du dv 



dudv 



+ t 



d 2 
\-f 

dt^ 



d__d_ 

du dv 



92 



9^ 9^ ^ 



dudv 
log: 



log I 



d_ d_ 

du dv 



logP 



(46) 



3 I a 

Ou dv 



^in terms of the Wronskian, defined in subsection 11.31 
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In the TO = 2 case, if U^^^ = {—oo,u), U^^^ = {oo,v), U^'^^ = {oo,w), ti = t and t2 = s the result for 
logF^^y {t,s,u,v,w) is 



Corollary 2. The logarithm of the 3-time joint probability for the Airy process F^"^{t, s, u, v.w) ( abbreviated 

du ' dv ' du-' 



as logPj satisfies a third order non-linear PDE in variables u, v, w, t and t (D = -§- + -§- + -§-) 



_ d_ d_ d_ 

du dv dw 



d d 
'dt+'dS 



+ t{s - t)D 
+{2t - s){2s - t){s + t) 



d_ 

dt 

g3 



£)logP + 



g3 



g3 



Q3 



dv?dv du^dw dudv"^ 



Q3 



dudvdw dudw"^ 
1 
2 



-{s-tf 



dv'^dw 



+ {s-tf 



g3 



dvdw'^ 



log I 



d . . d d 
s— + {2t - s)— + s — 
ou av aw 



£)logP,£>MogP^ . (47) 

D 



1.3 Notational convenience 



Throughout this paper, parentheses (...) always include numbers and functions; brackets [. . . ] always include 
operators; braces {. . . } are always for Wronskians: {/, gj/j = gDf — fDg, where Z) is a differential operator. 
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2 The joint probability in the Dyson Brownian motion 

To get the PDE, we need to consider a generalized integral such that indices i and j can be any positive 
integers 

r.(tf\cS;af) = l/.../ V{X^'^)V{Xi-^) 

^ J JuWx-xU(^)'' 

rn rn m n 



1=0 



1=0 k=l 
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with C a normalization constant such that pDyson _ t„L, where the locus £ is defined as = 1, 2, . . . , to— 1, 



1,2,, 



, m, Cfc 



V:(0) _ _ 



,(0 



Am) 



1 - 



1 



1 - 
2cfc 



(49) 



1 - c- 



Remark 2. All and cj';!- are variables of 
the value 0. Therefore it is legitimate to consider 



„(0 



and all other coefficients 0. 



in latter part of the paper, though most of them assume 



9t ' 



iDyson 



Remark 3. Since we allow s^'^ to be —oo and 03^^ to be +00, the integral in 



values of i and c)'",' . However, if we assume t. 



may be divergent for general 
for i > 2 , = for max(i, j) > 1, and values of t^'-*, 
and cl'^l are near to the locus £, then the integral is convergent, and all algebraic operations in latter 
part of the paper can be taken in this restricted setting, so they are legitimate. 



Now we consider actions of Z?''^ on Tn- Since Z)'^^ acts on the integral domains of X^^\ . . . , An'', by the 
formula 

d d 



f{x)dx, 



(50) 
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we get 



1 2ro 

-V 


d 


1=1 




m 




1=0 








hi- 


/ 


m 





C/(o)"x---x!7("')' 



1=1 



y(A(°))y(A('")) 

m 

;=o k=i 

( y(A(°))y(A(")) 



;=o 1=1 



m n 

(!) v-n \('-i)S(')M -n- -n- ,, (i) 

fe 

=0 fe=l 



1 



oo n 



oo 



n 



i=l fc=l id=0 k=l 



/ m 771 



\ T-r T-r (I) 



=0 /s=l 



1 

'"c 



(7(0) "x-'-x 



,(0) + ^^(1) + ^'^i.i ^^(1) 



i=2 



i=l 



„(0 xC-i)' xO' 



=2 j=l 



(0 



=0 k=l 



CO J5 ^ Q 

™1 +Z^"i „,(0) + 2^^1,i 



oc oc 



o (0) /l^ 

i=2 OT^-l i=l i=2 j=l '-''^i-1,. 

and similarly = 1, . . . , m — 1) 



r\ r\ OO OO 



nt 



dtti 



i=l 



EE^-^: 

j=2 i=l 



<=i 



5 



at. 

nt^+E^t 



(i+i) 



EE-S 



i=2 



at: 



i=2 j=l 

oc 

+ E^a 



(m) i9 \ - (m) d \ - \ - (m) 5 



(m) 



-1 i=l 



\ ^ \ ^ . (m 
o (m-1) +Z^Z^.^'=i,J „ (m) 

j=2 1=1 <^S,j-l 



On the locus £ we get (We abbreviate P^y^°'^ as P„ here and latter) 



2 a _^ 2ci 9 



(1) 



Now we define anm+lxm+1 matrix 



« eg 



J = 



„(i) 2/(1) 

^1,1 "^''2 ^-1,1 



'-1,1 

(m) Q,(m) . 
\ Cl,l ^^^2 / 

whose rows and columns are indexed from to m. On the locus 

2ci 



2ci 



l-d 



2ci 



1-c? 



i + l 2C2 



2C2 



1-c? 



2c„ 



2c„ 
2 



\ 1-3^ ~ 1-c^ 

We can find the entries of the first and the last row of J on the locus exphcitly. 

I 



Jm,l S, 



m—l 



1 



2 XX Cm-i+l — 2^ 
j=l 



tl—t-n 



and especially 



1 1 



i=l 



Then let 



y^im.iy 



/dO'1\ 

^1,1 

y£,m,iy 



we have 
Lemma 1. 



_g;0,l_g;m,l j p ^ £;m,i£;u,i j p ^ 



m.l 77'0,1 1 



9^ log I 



(57) 



(58) 



(59) 
(60) 



(61) 



(62) 



(63) 



Proof. First, since E°'^ and iS^'i are linear combinations of i^^'i's, they are differential operators of order 
1, and we have 

RO,lp Rm,lp__ RO,lpm,lp 

(64) 



^0,l^m,l ^Qgp^ ^ n_ n _^ 



P2 
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By dSl - (EH), & and §UI, we get 



(65) 



and 



Therefore 



= J2 J-^.^D'^^ = 



\i=0 



at'"' 9t( 



(m) ri 



(66) 



(67) 



Here we need to be careful about the term — j^^n- By ([5T|) - ([55]) . the action of £'™'i on t„ is 



ati 



equivalent to that of a differential operator which does not contain o'p explicitly. On the locus £, all terms 
of the differential operator except for ^^^^ vanish, so we can ignore them and replace E"^'^ by ^^f^^ between 
E^'^ and P„. 

Since E°'^ and commute, 



ati 



0,1 d 



d 



:E" 



by (in]) - ([531), (ISll) and ((621), we have the identity for the action of E^-^ on r„ 



aii" 



(68) 



E^'^Tn =Jo,0 



nt 



(1) 



d 



i=2 i=l 



dt 



(1) ' EE*4i o (1) 



(=1 



n4'^ + E 
(9 



9 



oo oc 

+EE< 



dt 



E^S 



i=2 ] = 1 

d 



OO oo 



5 



dt 



j=2 i=l "'^1,1-1 



d 



i=2 j = l 



(m) 



(m-1) 



EE 



(m) 



j=2 1=1 '^'-ij-1 



(69) 



with coefficients of all terms except for ^^fm) + "■■/o.mii™'' of the right-hand side operator vanishing on the 
locus £ and not containing 4™^ explicitly. So on the locus 



d 



?0M 



(m) " 



? — 7 I 7 — T + TT-Jq mt\ ^ + 



and 



i;"'i£"'iiogP„ = 



ap„ ap„ a^p„ 
at^"' atj"' , at 



(70) 

(71) 
□ 
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Similarly to ([?T|) - with the help of the formula 



we get ( 



2^k=l g^(0} Afe 



, 2ro 



is regarded as an operator) 
(0) d 



da 



(0) 



(7(0) "X---X [/("•)' 



y(A("))y(A('")) 



1 



1=1 



1=0 k=l 



C.J ./(7(0)"x---x(7(") 



^ 9 AO) 

k 

(0 1,(0* -I-r ^(0 



.fe=l CAj. J \ 



1=0 1=1 



( F(A(''))y(A('")) 

m n 

UUdx 



1 

"c 



oo n 



ij=0 /c=l 



;=0 k=l 



i=l k=l 



/ rn rn 



m n 

i<'-i)*x<')J \ T-r T-r (/) 



1 

c 



1=0 



[/(")" x---x;7(" 



!=0 fe=l 



(1) d n{n + 1) 



^ (0) ^ »J ^ (1) 



/ m m 



m 71 

i(<-i)'x<')J \ T-r T-r (/) 



/=0 

C30 oc 



E^^f'T^^+EE 



f)f(0) ' ''3 ^ (1) 



1=1 

. (I) d n(n + 1) 



;=0 k=l 



and similarly 



E 



it 



(m) d 



dt 



(m) 



EE^- 

i=i j=i 



(m) d n{n + 1) 



^3 



and for I — 1 , . . . , to — 1 



5 



EE^^Stw + EE^ 



1=1 j=l 



= 1 3 = 1 



Jl+1)_0_ 
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On the locus £ we get = 1, . . . , m — 1) 



d 



2ci d n{n + 1) 



D 



l,2jl 



^l+l 



d , 2ci d , 2Q^ 



9 



2 9 2cm 9 n(ri + 1) 



If we define {I = 1, 2, . . . , m — 1) 

pi. 2 _n''2 J') ^ 

^m,2 2 _ ^(™) 1^ 



(76) 
(77) 
(78) 

(79) 
(80) 
(81) 



we have 



Lemma 2. For k,l — 0,1, ... ,m, 

E'^'^E'-'logP^ = 2ti%-^^ + 5iE^^HogF„. (82) 

Ct'2 (yi^-^ 

Proof. With arguments similar to those for (164^ and ([67|) . we get 



s P2 P2 p^^ V ; 

Here 

-^5 — - fc = 1 or TO, 
C={ ^ (84) 



— n 



otherwise. 



Similar to and we have 



and 



with coefficients of all terms other than ^J^) , ^^^T ^ right-hand side operator vanishing 
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on £ and not containing t^''' explicitly. Therefore with an argument similar to that for (j70p 

'24'=) + C¥„) ^ 24^) M d^c] 



since 



Since on the locus £, c"^l and 4''' are functions of ci = e . . . , Cm = e defined in (|49ll . by the 
rule we get as operators on P„ (/ = 1, 2, . . . , m) 

d _ 2cf d 2c2 a 2q(1 + c2) d 



and 

m a _ 2cf 5 ^ 2cf d 1-cf d 



'■''dc['\ ^-cfdtt'^ l~cf94) ^ + ^fdsi 
Therefore hy ^ - ^ we get on £ that (/ = 1, 2, . . . , m - 1) 



^0,2 ^^0,2 



2c2 d 2cl d 1-4 d 



^m,2 j^ni,2 



2 

2ci 9 2ci 9 1 - d a 



Now we denote (/ = 1, 2, . . . , m — 1) 



^0,2 ^^0 



.2 



i-cl d 

1 + c? 9si ' 



jjm.2 



l + cf dsi 1 + cf_^_-^ dsi+i ' 
1-ci d 



1 + c2 9s„ 



and we have 
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Lemma 3. For ^ = 1, . . . , to — 1, 



F"'2£;™^ilogI 



92 



2cl 



52 



2c? 



92 / 2 



logP„, 



2ci 92 



2r4 



logP„ 
2 



92 



Wail™) 



^2 



Proof. We only prove (|97p . By (pij) . on £ we have 

2c^ 9 2c2 a 



log I 



^0,2 ^ ^0,2 



and similar to (p7)) and ([701) . we have 

9 



S''ilogP„ = - 



so that with the result of 



(fc) 



dt). 



F°'2£;'"'ilogP„ ^E"'^E"'^HogF,, + -M__|_£;™4 



=2t 



(0) 92iogI 



1 - dt^ 
2cl 92 log J 



log I 



9'logPn 



1 - 4 941) 



2c? 92 log J 



2c2 



logP„. 



Finally we define 



where 



K 



l-cf 



/G0-2\ 
^1,2 

yG"^2y 



2c^ 



/f0-2\ 
^1,2 



2ck 2ci 



1-4 1-4 1-4 
2c| 



24 



2 



15 



We can get K \yy substituting each q in J ^j^ by cf, so we have 



1=1 



1 ' 



m — l 



Kmd = "2 lie: 

i=l 



In 



2 

m—i+l ■ 



and get by (HH) - ([99 
Lemma 4. 



G°'2i;"'ilogI 



logP„ + Xo,™i^"'MogI 



Symmetrically, we can get by the same method 
Lemma 5. 

G"^2£;^^MogP„ = 



92 



logP„ + /f™.o-B°''log 



By the result of P, 



5 . r„+i 
log ■ — 



(0) 



d 

Am) 



Tn-1 



92 



— log r„ ' 



log 



T , -, aj(»")aj(0) l0gT„ 



we get the differential equation on the locus 



0^ 



log T„ ' 



Tn-l 



E""^^ log ■ 



Tn-1 



G°-'-E"'^l0gT„|^- j^O.n.E^-llogT^I^ 

£~ £;o,i£;m,iiogr„|£ -nJ„,,o 

_G"'2^°'ilogT,4 - j^„,o^"'ilogT„|^ 
^- ^™'i£;o.ilogT„|£-nJo,„. 



where the result 



can be proved like lemma [TJ By the identity 



^0,l^m,l log !»±i 



i^o-ilog 



Tn+1 



Tn-1 



Tn-l 



= s™'ii;"'iiog^ 



Tn-l 



we get the final result 



^0^1 G^'^gO'l l0gf„ - g,n,O^P-l l0gf„ ^ ^„ i G0.2^"'llogf„ - i^o,r„£^™''logf„ 



£;,ri,i£:o,ilogf„ -nJo,, 



SO,l^m,llogf„-nJ„,,0 



Now we denote 



Ai = 






Si = 






^2 = 


- 2(G"'2 - 




62 = 


- 2(G"''2 





(105) 



(106) 



(107) 

(108) 
(109) 

(110) 
(111) 

(112) 
(113) 
(114) 

(115) 
(116) 
(117) 
(118) 



and we get the equation ([34] 
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Remark 4. In the 2-time case, i.e., m = 1, 

Ai^D°^^ +ciD^'\ (119) 
Bi ^ciD°'^ + D^'\ (120) 

A, + clF'-' = i^O'i + clD'-' + (1 - c?) A _ c?, (121) 

Oti 

B2 ^cjF"^^ + Fi'i = cjD°-' + D^-^ + (1 - cl)-^ ~ c\. (122) 
Our PDE (1211) agrees with that in 12]. 

3 The joint probability in the Airy process 

In this section we adapt notations defined in (P^ll - pSl) . and by remark ITJ a'^ = a'^ = — oo. We denote 

(/=0,l,...,77l) 

2ri „ 2ri „ 

if all c^r} < +00, otherwise drop the part, we can write the diff'erential operators defined for the Dyson 
process as 

m 

y^i =v^^e-*'/"£'''\ (124) 

(=0 
m 

6i =v^^e(*'-*'")/"L)''\ (125) 

1=0 

m m m 

A2=Y^ e~2tVn^L2 ^ 2fi2 ^ e-2t,/"^M + n^(l - e""*'/")^ - e-^'-/^, (126) 

1=0 1=0 1=1 ' 

m m m 

^2 =^e2(''-*'")/"5'^2 _^ 2n2^e2(f,-t„)/n^i,l ^ -^(g2(t,^f„)/n _ ^~2t^/n)^ _ g-2t„/f.^ (^27) 

/=o ;=o 1=1 ' 

It is not difficult to see that (IMI implies 

[^ii32^i - 61^2^1] logPn • (AiBi logP„ + 2ne-*'") = 

Sa^i logP„ • AiBiAi logP„ - A^i logPn • BiAiBi logP„. (128) 
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Take substitutions P^ - P?7|) into we get 



L/=0 



i=0 



i=0 



9 



1=1 "^^ 



.1=0 



,^(l_e-2*Vn)^_e-2*™/" 



1=1 



d_ 

'dil 



n 

logP„ 



logP„ 



i=0 



/=0 

m 



.i=0 



^f,{ii-ir„)/njjl,l 



1=0 



logPn + n'^e 



2(t,-t„)/fi _ g-2t„/fi 



1=0 



1 = 1 



^g-t,/n^i,l 



.1=0 



logP„ X 



,/=0 



.1=0 



logP„ 



d 



J2 e-2tVn^/.2 ^ 2n2 ^ e-2*V"^M + n^(l - e'^''/")^ - e 



-2tm/fi 



L/=0 



.i=0 



logP„ X 



^e(ti-t™)/"5',i 



J=0 



^g(t,-*™)/nj3)i4 



.;=o 



log I 



Since we have commutator formulas 



n{e 



2{ti-im)/r, 



1=0 



d 



i=l ' 1=0 



m 
1=0 



(ti-2t„)/n _ ^-{ti+2t^)/n^j=)l,l 



(ti-3t,„)/n _ ^-(t,+t™)/n-j^!,l 



.Z=0 



/=0 



^g(t;-t,„)/«2)'.\ ^g2t,/n^i; 



.Z=0 



/=0 



= ^g(tl-2t™)/ft2)'.\ 
m 

^^g-(t!+t™)/n2)i,l^ 
i=0 



(129) 



(130) 
(131) 
(132) 
(133) 
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we can write (|129p as 



H 2 



L;=o 



,2(t,-t,„)/fi 



J=0 
H 2 



1=1 



^g(ti-*™)/fif)',i 
Li=0 



.i=0 



!=0 



'dt. 



logP„ 



logP„ 



;=o 

m 

^g(t,-2t,„)/n^;,l 



J LZ=0 



logP„ + n^e 



2p-t™/fi 



i=0 



i=0 



log I 



,i=0 



logP„ X 



.1=0 



^e(t!-*™)/ft£)',i 



log I 



U=0 



e-2tV"5M _^ J2 e-2*Vft5;,i + n^(l 



^g-(t(+t™)/n2)M 
ll=0 



J Li=o 
logP„ I X 



1=0 



1=1 



'dL 



logP„ 



Li=0 



H 2 



J li=o 



log I 



(134) 



Since all terms of the PDE involves n, we can expand the PDE with respect to n, with formulas (* can 
be 1 or 2) 



m m - m ^ m 



z=o 



1=0 



i=0 



m m c)"^ i-)^ /l^^ /1\ 

■m m ^ m ?n 



(135) 
(136) 



1=0 



1=0 



1=0 



1=0 



(137) 



I I L lit. ^ 1 1 L ^ I 11 

1=0 1=0 1=0 1=0 



4 



1=0 ^ ^ 



(138) 



1=1 



d_ 



2E^'^ + -EMt.-2U^ + 0(^^j, 

/—I /—I '' \ y 



1^1 



1 = 1 1 = 1 

d 2 9 



' dti n ^ dti 

1=1 1=1 



(139) 
(140) 
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Although the left hand side of (|129p contains 0{n'^) terms, after careftil calculation we find all 0{n'^), 0{n^) 
and 0{n^) terms disappear, and the equation becomes 

n 2 



1=0 



5: DM 

D 

m 

m 

.1^0 
m 



.1=0 



1=0 



1=1 



dti 



Y{2ti - 



d_ 

dti 



.1=0 



.1=0 



.1=0 



logP„ 
log I 

logE 



.1=0 



1=0 



^(2f, - tra)&'' 



.1=0 



logP„ 



.1=0 



Y&^^ 



.1=0 



logP„, 



Y&^^ 



.1=0 J 11=0 



Y&^' 



logP„ 



n 



The term O (i) in plT|) is a quadratic function in term of logP„ and its derivatives with coefficients 
O (4). By the definition of F-^"y in (pS)) and the convergence result in we take the limit n ^ oo, and 
get the PDE (|45|) after the changing of notations, i.e., cleaning all "bars" for variables and operators. 
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